Soliton Propagation in Chains with Simple Nonlocal Defects 
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We study the propagation of solitons on complex chains buih by inserting finite graphs at two 
sites of an unbranched chain. We compare numerical findings with the results of an analytical linear 
approximation scheme describing the interaction of large-fast solitons with non-local topological 
defects on a chain. We show that the transmission properties of the solitons strongly depend on the 
structure of the inserted graph, giving a tool to control the soliton propagation through the choice 
of pertinent graphs to be attached to the chain. 



I. INTRODUCTION 

Recently, much attention has been devoted to the anal- 
ysis of complexity arising in discrete physical systems 
living on networks with non-trivial topologies: remark- 
able examples are by now given by networks of nonlin- 
ear waveguides [1], Bose-Einstein condensates in optical 
lattices [2], Josephson junction networks (JJN) [3] and 
silicon-based photonic crystals [4]. In these systems, the 
choice of the network's topology allows to engineer new 
macroscopically coherent quantum states: a remarkable 
example of macroscopic quantum coherence topologically 
induced by the network topology in JJN has been pre- 
dicted in [5]. 

Another relevant area where one should be able to ev- 
idence new phenomena induced in discrete quantum sys- 
tems by the topology is provided by nonlinear dynamical 
systems. Nonlinearity already produces remarkable phe- 
nomena such as soliton propagation [6] even in very sim- 
ple geometries. A few steps in the study of the interplay 
between nonlinearity and complex topology have been 
made: recently, the effects of uniformity break on soli- 
ton propagation [7,8] and localized modes [9] have been 
investigated by considering y-junctions [7,8] (consisting 
of a long chain inserted on a site of a chain yielding a 
star-like geometry) or geometries like junctions of two 
infinite waveguides or the waveguide coupler [9]. In [10] 
it has been considered the discrete nonlinear Schrodingcr 
equation (DNLSE) on a discrete network obtained by 
inserting a finite graph at a site of a one dimensional un- 
branched chain and the soliton propagation through this 
finite graph has been studied by numerical and analyti- 
cal tools. It has been showed that for sufficiently large 
and fast solitons the soliton momenta for perfect reflec- 
tion and transmission can be analytically related to the 
energy levels of the inserted graph. Such results for the 
transmission properties of solitons in inhomogeneous net- 
works have been used in [11] to show that it is possible 
to engineer topological filters for the soliton motion on a 
complex chain. 

The simple criterion relating the perfect reflection and 
transmission momenta and the energy levels of the in- 
serted graph obtained in [10] has been obtained in the 
situation in which the graph is inserted at a single site. 



In the present paper we address the more difficult is- 
sue of the soliton propagation for the DNLSE on a net- 
work obtained by inserting a graph on two sites of an 
unbranched chain. In this respect, the inserted graph 
has an internal structure and it is seen by the soliton like 
an extended topological defect, since during the motion 
along the unbranched chain the soliton wavefunction is 
modified in more than one site. One therefore expects a 
variety of possible resonances between the energies char- 
acterizing the soliton propagation and the energy levels 
of the network. We shall focus here only on some partic- 
ularly simple inserted graphs: loops, bubbles and single 
links attached in two sites (see Fig.l). The plan of the 
paper is the following: after introducing in the next Sec- 
tion the DNLSE on a graph, we report in Section III 
our results for the transmission properties of solitons on 
such networks by resorting to numerical simulations and 
analytical results based on the linear approximation for 
the analysis of the interaction of fast solitons with these 
topological defects. Section IV is devoted to our conclu- 



II. THE DNLSE ON A GRAPH 

The DNLSE is a paradigmatic example of a nonlin- 
ear wave equation extensively studied on regular lattices 
[12-14]. On a chain it reads 
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(1) 



where n = • • • , — 1, 0, 1 • • • is an integer index denoting the 
site position and A is the coefficient of the nonlinear term. 
The normalization condition is | ■(/'« 1- It is well 
known that the DNLSE on a homogeneous chain is not 
integrable; however, soliton-like wave-packets can prop- 
agate with (quasi)momentum k for a long time [15]. By 
means of a variational approach for gaussian wavepack- 
ets with width 7 much larger than 1 (the distances are 
in units of the lattice length), one finds that, for A > 0, 
it is possible to have solitonic solutions of the variational 
equations of motion with constant width 7 if cos fc < 
and A equals the critical value [16] 
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Numerical simulations confirm that the stability of these 
wave packet is robust for long times. 

The generalization of Eq.(l) on an arbitrary graph is 

^^ = -^EAi^^^+^IV'd'^.: (3) 

j 

in Eq.(3) Ai^j is the so-called adjacency matrix of the 
graph [17], which is defined to be 1 if i and j are nearest- 
neighbours sites, and otherwise. We shall limit ourself 
to networks obtained inserting simple graphs at two sites 
of the unbranched chain (Fig.l). The two sites of the 
unbranched chain at which the graph is inserted are de- 
fined to be n = and n = n. We assume that the soliton 
is traveling from the left with constant velocity v, re- 
lated to fc by w ~ sin k. Eq.(3) is numerically solved with 
A = Ksoi using as initial condition 

= 0) = /Ce-("-«'')'/^'+^'=("-^°) (4) 

at the sites n = -- -,— l,0,l---of the unbranched chain 
and = 0) = at the sites of the added graph. In 
Eq.(4) /C is a normalization factor and is the initial 
position of the soliton center: we choose < with 
I ^0 1^ 1 and 7r/2 < fc < TT, so that t; > and the soliton 
moves from the left to the right of the unbranched chain. 
From the numerical solution at very large times (well af- 
ter the collision with the inserted graph) the reflection 
and transmission coefficients R and T are compiited by 
the relations R = X)„<o I P and T = Y.n>n I "^n 1^- 
In the following we shall present numerical results for 
the coefficients R and T obtained (for an initial width 
7 = 40) for different values of k and for the networks 
plotted in Fig.l. 

When the soliton is large (7^1) and fast enough that 
the soliton-defect collision time is much shorter than the 
soliton dispersion time (i.e. the time scale in which the 
wavepacket will spread in absence of interaction), thus 
one may resort to a linear approximation to compute 
the transmission coefficients [18,19] since, in these limits, 
the soliton may be considered as a set of non interacting 
plane waves experiencing scattering on the graph. The 
soliton transmission coefficients may be then estimated 
by computing in the linear regime the transmission co- 
efficients of a plane wave across this topological defect 
inserted at two sites. Afterwords we shall compare our 
analytical findings with the results coming from the nu- 
merical solution of Eq.(3). 

III. TRANSMISSION COEFFICIENTS 
A. Loops 

Let consider the situation in which a loop with length 
L is inserted at the two neighbouring sites n = and 



n = 1 of the unbranched chain [see Figs.l(a,b)]. In the 
linear approximation, which is expected to be reasonable 
when L < 7, the transmission coefficients may be deter- 
mined by considering a plane wave solution having 

V)„ = ae''="-F6e-*'=" (5) 

for n < 0, while for n > 1 one puts 

Vn=ce^'=". (6) 

At the sites a = 1, • • • , L of the loop the wavefunction is 
given by 

=rfe''="-F/e-''=". (7) 
The eigenvalue equation to solve is 

where i and j run on all the sites of the whole network. Of 

course, from Eqs.(5)-(6) one obtains £ = — cosfc. From 
the continuity in n = 0, n = 1, a = 1 and a = L one 
gets, respectively 

- i (a e-'^ + be'^ + c e'*^ + de'^ + f e''^) =£{a + h) 

(9) 

-i(a + 6 + ce2*'= + (ie''=^+/e-*'=^) =fce*'= (10) 

a + b = d + e (11) 

ce'^ = d e'^'^^+i) + f e-^'^^^+i) (12) 

Prom Eqs.(9)-(12) the reflection coefficient R =\ b/a p is 
given by: 

^ _|_ p2tfeL _ p2ife(L + l) _ 2gifc(-t+2) 
^ ~ I _ 3g2ifc _j_ g4i/c _ 2e»'=(-^^+2) + g2ife(L+2) _|_ 2e^HL+4:) 

(13) 

i.e., R = Rm /Rd; where Rm = 5 — 4cos(2fc)-|-4cos(fcL) — 
cos(2fc(L+l))-4 cos(fc(L+2)) and Rd = 10[l-cos(2A:) + 
cos(A;L) - cos(fc(L + 2))] + cos(4fc) - 2cos(A;(L - 2)) + 
cos(2A:L)-3 cos(2A;(L+l))+cos(A;(L-t-2))+2 cos(fc(L+4)). 
Of course, if one inserts a loop at a single site of the un- 
branched chain (i.e., n = n), one has to require 'ipa=i = 
'^a=L and Eqs.(ll)-(12) simply become a + b = c = d-\- f . 

The results obtained for the reflection coefficient R 
from the numerical solution of the DNLSE (3) and from 
the linear approximation, Eq.(13), are reported in Fig. 2 
for L = 2 (circles) and L = 4 (crosses) and are in re- 
markable agreement. We see that for i = 4 one has 
two values of the momenta for which there is perfect re- 
flection [k « 2.10 and k « 2.25): this shows that one 
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can control the soliton propagation by properly choos- 
ing the topology of the network. The average position 
(") = Sn'^ I V'n(i) P is plotted vs. time in Fig. 3. 
We see that before and after the collision the soliton 
move with constant velocity: with v = d{n)/dt, one has 
V ~ sin(fc) before the collision, and v ~ — sin(A:) after the 
collision which (almost) totally reflects the soliton. The 
splitting of the soliton in transmitted and reflected parts 
is illustrated in Figs. 4-5, where we consider a loop with 
length L = 2 and k = 1.8. In Fig.4 we plot | tpn P 
at five different times, including a time {t = 200) in 
which the soliton hits the loop. In Fig. 5 wc plot the 
time evolution of the number iV; of particles in the left 
{^i = Z^n<o I V'n 1^)) the number of particles in the 
right {Nr = J2n>i I V'n 1^) ^ud the number Nioop of par- 
ticles in the loop {Nioop = J2a=i I ^« 1^)- ^^^^ ^^^^ 
for t w 200 the number of particles in the loop increases 
and after th(^ reflection decreases. 

In the simplest case (L = 1) a single extra site is at- 
tached to the sites n = and n = 1 of the unbranched 
chain, as in Fig. 1(a): then Eq.(13) simplifies to 

2cos^(fc/2) 
2-|-cosA;-cos(3A;)' ^ ^ 

The comparison between the results for R from the nu- 
merical solution of the DNLSE (3) and from Eq.(14) is 
reported in Fig. 6, showing also in this case a good agree- 
ment. 

We observe that the extra sites of the inserted graph 
can be view as external Fano degrees of freedom coupled 
to the chain [19-21]. In particular, in [21] an additional 
discrete state is coupled to the sites of a straight linear 
chain: this would correspond in our description to a site 
linked to all the sites of the unbranched chain. 



B. Bubbles 

We refer in this Section to inserted p-bubblc graphs, i.e. 
to loops with length L = 1 inserted j3-times at two sites of 
the unbranched chains which are distant 2 [see Fig. 1(c)]. 
To fix the notations the sites of the unbranched chain are 
defined to be n = • • • , —1, Oi, 1, • • •, and other p — 1 sites 
O2, 03, • • • , Op are linked to the sites —1 and 1. Of course, 
if p = 1, we have the simple chain. 

The coefficient R may be obtained in the linear ap- 
proximation in the following way: we assume ipn = 
ae*'^" + 6 6"*'=" for n < -1, Vn = ce**=" for n > 1, 
and tpo = fpOi = • • • = (for symmetry, all the 

sites inside the bubble are equivalent). The eigenvalue 
equation (8), with £ = —cos A;, in the site —1 reads 
-(l/2)(^_2+pVo)=<?V'-i,i-e., 

-\{ae-'^'^ +he^''' +pi)Q) = £{ae-'^ +he'^). (15) 

Similarly, Eq.(8) in the site -|-1 gives 



-\{p^Q + ce'''^)=£ce-'^. (16) 
In a site inside the bubble, Eq. (8) reads 

(oe-*'=-h6e*'=-hce-*'=) =£:^o- (17) 
Prom Eqs.(15)-(17) one gets 

R= ^^72 (18) 

I-h(^) tan2/c 

of course, when p = 1. no reflection occurs. Eq.(18) 
shows that increasing k from 7r/2 to tt (i.e., decreasing 
the velocity), the reflection increases: slower solitons are 
more reflected. In this meaning, the p-bubble is an high- 
pass, i.e. only solitons with high velocity are transmitted: 
by varying p one can control the width of the range of 
transmitted velocities. In Fig. 7 we plot the reflection co- 
efficient R for p = 2 and p = 20, showing that a larger p 
make smaller the range of transmitted velocities. As in 
Figs. 2 and 6, solid lines correspond to the analytical esti- 
mate [given by Eq.(18)], and numerical results from the 
DNLSE (3) are expressed by circles (p = 2) and crosses 
(p = 20). 

C. SepEirate links 

In this Section we consider the effect on the soliton 
propagation of two extra sites linked to two sites the 
unbranched chain, as in Fig. 1(d). To fix notations, we 
suppose that the two extra sites, a and /3, are linked 
respectively to n = and h = L. 

To obtain the coefficient R in the linear approximation, 
one proceeds as before: one assumes V'n = a e'*^" -|- b e"*'^" 
for n < 0, V'n = c e**^" for n > _L, and ip^ = d e^^" + f e**^" 
for n = 1, • • • , L — 1. The eigenvalue equation (8), with 
£ = — cos k, in the sites 0, L, A, B, 1 and L + 1 reads 
respectively 

~ {ae-'^ + be'^ +i)oc+de'^ + fe-'^) =£{a + b) 

(19) 



(20) 

-i(a + 6)=£:Va (21) 
-^ce"'^=£^f, (22) 
a + b = d + e (23) 
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ce 



ikL ^ ^^ikL _^ j^^-ikL_ 



(24) 



One has seven unknowns (a, 6, c, d, /, tpa said ijjp) and 
the six equations (19)-(24) (the remaining is provided by 
the normalization). Solving for b/a one gets 



R = 



1 - e^*'^ - e 



„2ifcL 



, g2jfc(L+l) _ g2ife(L+2) 



g2ife(i+2) _ ^g2ifc _|_ g4ife _ J^'^^ 



for L = 1 Eq.(25) simplifies to 



R 



(l + 2cos(2fc))^ 



9 cos2 k + (sin fc + 2 sin (3A:)) 
and for L = 2 to 



R 



(1 - 2cos(2A:)) 



4sin^ (2fc) + (2cos(2fc) - 1) 



2 • 



(25) 



(26) 



(27) 



Notice that if you have only a site linked to a site of the 
unbranched chain, one obtains in the linear limit [19,10] 



R 



1 



l + 4sin^(2fc) 



(28) 



a comparison of Eq.(28) with Eq.(25) shows the remark- 
able effect (also for large L) of the second added link on 
the transmission properties. In Fig. 8 we plot the reflec- 
tion coefficient i? for L = 1 and L = 2, showing that 
in both situations one has a peak in the transmission 
{R w 0) for two different values: the two separate links 
behave approximately as transmission filters, i.e., allow- 
ing for the transmission of solitons with certain velocities. 
Increasing the distance L between the sites at which the 
added sites arc linked the number of such transmission 
peaks in turn increases. In Fig. 8 the solid lines corre- 
spond to the analytical estimate [given by Eqs.(26)-(27)], 
and numerical results from the DNLSE (3) are expressed 
by circles (L = 1) and crosses (L = 2). 



IV. CONCLUSIONS 

In this paper we studied the propagation of large-fast 
solitons on complex chains built by inserting finite graphs 
at two sites of an unbranched chain. The inserted graph 
is seen by the soliton like an extended topological de- 
fect: the transmission properties of the solitons strongly 
depend on the structure of the inserted graph. We con- 
sidered simple inserted graphs: loops, bubbles and single 
links. For inserted loops peaks for perfect reflection oc- 
curs, while for bubbles the inserted graph behaves as a 
high-pass filter. When two added sites are linked to two 
sites having distance L between them, transmission peaks 
appears and increasing the L the number of such trans- 
mission peaks increases. In all these situations we com- 
pared numerical findings with the results of an analytical 
linear approximation, obtaining a good agreement. 



In conclusion we think that the study of nonlinear dy- 
namical systems on complex networks is a wide subject to 
investigate, and that the main motivation of such study 
is that the network topology provides a natural tool to 
control the nonlinear dynamics of wavepackets. In par- 
ticular, we mention as possible interesting future studies 
the study of the propagation of very localized breathers 
on complex networks and the nonlinear trapping of soli- 
tonic solution in chains with topological defects. 
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FIG. 1. Inserting simple graphs at two sites of a linear 
chain; (a) a single extra site attached to two sites of the 
chain; (b) a loop with length L; (c) a bubble with p = 3 sites; 
(d) two single links attached at two sites with distance L = 3. 



FIG. 3. Center of mass position (n) vs. time for k = 2.1 

and a loop having L = A, which corresponds to perfect reflec- 
tion. The solid line correspond to the numerical solution of 
Eq.(3), and the dashed lines to free motion of the soliton with 
absolute value of the velocity \v\= sin/c. 
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FIG. 2. Reflection coefficient i? as a function of k (with 
k between 7r/2 and tt) when loops with length 2 and 4 are 
attached. Empty circles {L = 2) and crosses {L = 4) corre- 
spond to the numerical solution of Eq.(3): as initial condition 
we choose a Gaussian with initial width 7 = 40 and momen- 
tum k (see text). Solid lines correspond to the analytical 
prediction (13). 



FIG. 4. Soliton propagation obtained from Eq.(3) [for 
momentum fc=1.8 and width 70 = 40] through an in- 
serted loop of length 2. The soliton profile is plotted for 
z = 0, 100, 200, 300, 400 corresponding to (1) • • • (5). 
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FIG. 5. Time evolution from the numerical solution of 
Eq.(3) for the number of particles on the left of the loop (Ni), 
on the right of the loop (Nr) and on the loop (Nioop) for 
L = 2 and k = 1.8. Around t = 200 particles enter the loop 
(compare with the previous figure). 



k 

FIG. 7. Reflection coefficient ii as a function of k when 
p-bubbles with p = 2 and p = 20 are attached. Empty cir- 
cles {p = 2) and crosses (p = 20) correspond to the numerical 
solution of Eq.(3), solid lines correspond to the analytical pre- 
diction (18). 




FIG. 6. Reflection coefficient _R as a function of k (with k 
between 7r/2 and n) for a loop made of a single site (L = 1). 
Empty circles correspond to the numerical solution of Eq.(3) 
and the solid lines correspond to the analytical prediction 
(14). 



FIG. 8. Reflection coefficient _R as a function of k when 
separate sites are linked to sites distant L = 1 and L = 2. 
Empty circles {L = 1) and crosses {L = 2) correspond to 
the numerical solution of Eq.(3), solid lines correspond to the 
analytical predictions (26)- (27). 
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